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In this paper the accelerating expansion of our universe at the late cosmic evolution time in a 
generally modified (extended) Chaplygin gas (Dark Fluid) model is detailed , which is characterized 
by two parameters (m, a). Different choices to the parameters m and a divide this model into two 
main kinds of situation by different properties with cosmological interests. With proper choices of 
parameters, we find that this extended model can realize the phantom divide w = — 1 (Equation 
Of State parameter) crossing phenomenon with interesting ranges of the scale factor value a{t) 
corresponding to w = — 1 and present value of state parameter w. 

Additionally, through Taylor series expansion of the function A(a) in the extended Chaplygin gas 
model , a specific Equation of State is gained. Under the framework of Friedman- Robertson- Walker 
cosmic model, it can successfully explain the accelerating expansion of our universe. However, the 
value of w in this case is large than —1, that is, indicating it like a quintessence fluid and no w — — 1 
crossing occurs. 

PACS numbers: 98.80.Cq, 98.80.-k 



I. INTRODUCTION 

Dark side of the Universe has been puzzling us across 
the century pj, especially the recent years discovery as 
coined Dark Energy. Observations of type la super- 
nova(SNe la) directly suggest that the expansion of the 
universe is accelerating, and the measurement of the cos- 
mic microwave background (CMB) and the galaxy 
power spectrum for large scale structure Q indicate that 
in spatially flat isotropic universe, about two-thirds of 
the critical energy density seems to be stored in a so 
called dark energy component with enough negative pres- 
sure responsible for the currently cosmic accelerating 
expansion^. It is clear from observations that most of 
the matter in the Universe is in a dark (non-baryonic) 
form (see, for instance, |5j). To understand the mysteri- 
ous dark components is even great challenge for tempo- 
rary physicists from this new century crossings. 

The simplest candidate for the dark energy is a cos- 
mological constant A, which has a specially simple pres- 
sure expression pa = — PA- However, the A-term requires 
that the vacuum energy density be fine tuned to have 
the observed very tiny value, the famous "old" cosmo- 
logical constant problem. So, many other different forms 
of dynamically changing dark energy models have been 
proposed instead of the only cosmological constant in- 
corporated model. Usually, the equation of statc(EOS) 
for describing dark energy can be assumedly factorized 
into the form of pde = wpde, where w may depend 
on cosmological redshift z [8j, or scale factor a(t) and so 
on. By the way, the case for w = — 1 corresponding to 
the cosmological constant, which has involved singularity 
in perturbation calculations, was thought as border-case 
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called the phantom divide 

To further explore the properties of the so called dark 
energy or its concrete form of EOS w, some authors pro- 
pose kinds of models such as px = —^-(Chaplygin gas 
model|I3), which is not consistent with astrophysics ob- 
servations as it can produce oscillations and exponen- 
tial blow-up to the matter power spectrum, its general- 
izations [Tfj like px = — ^-(generalized Chaplygin gas 

model(GCGM)@,Jl2), or generalizing the constant A 
to variable forms[6|, Px = ce(px — Po) (linearized EOS 
model0|, compared with the conventionally perfect fluid 
EOS px = (7 — an d so on(subscript X represents 

dark energy fluid except for the cosmological constant). 
Among these models, the fluid of GCGM has a dual be- 
havior: it mimics matter(p ss 0) domination at very early 
stage in the history of the universe evolution and a cosmo- 
logical constant much later(with a smooth transition in 
between) , which is highly suggestive of a unified descrip- 
tion of dark matter and dark energy (see reference 
that we may address it as Dark Fluid, as reason above. 
Moreover, dark energy may be made up of different en- 
ergy fluids with different EOS, and these different energy 
fluids interact with each other through p' s — Ps p ^ p p 3 = ^- , 
where subscript s indicates a component , which is called 
the self-interacting fluids(see reference \u\)- 

With the merits of GCGM, we further discuss a mod- 
ified general Chaplygin gas model by extending its EOS 
as an explicitly scale factor related form, which is based 
on the original relation px = Moreover, through 

the mathematically Taylor series expansion of the rela- 
tively general EOS we elaborately demonstrate concrete 
values of two characteristic parameters (m, a) for this 
ECG model, and we find the cosmic effective composi- 
tions (matter, radiation and vacuum energy) are from 
the extended Chaplygin gas component. 

The plot of this paper is arranged as: In Sect[H] first we 
introduce the framework by mathematically integrable 
expression of the extended Chaplygin gas model with an- 
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alytical investigations. From the mathematical point of 
view, we divide the model into two classes described by 
two characters m and a. And then in Subsects.IIA and 
IIB, we investigate the two classes with general theoreti- 
cal discussions, respectively. Secondly, through the Tay- 
lor series expansion method , a set of parameters (m, a) 
is determined and its physical indications are presented 
in Sect IIIII At the last part, we summarize with discus- 
sions. 



II. EXTENDED CHAPLYGIN GAS MODEL 

The metric of a homogeneous and isotropic universe is 
usually written as follows 



ds 2 



dt 2 - a 2 {t)dl 2 



(1) 



where dl 2 is the metric of a 3-manifold of constant cur- 
vature (k = 0, ±1), and a(t) is the expansion factor. 

Under the framework of Friedman-Robertson- Walker 
cosmology, the dynamic evolution of cosmology is com- 
pletely determined by the Friedman equation 



k 



, energy conservation equation 

d(pa 3 ) = -pd(a 3 ) 



(2) 



(3) 



and equation of state p = f(p) if we have. With the con- 
sideration that our Universe is filled mainly with both 
matter (mostly non-relative dark matter) and dark en- 
ergy( the radiation component negligible at the later 
phase for the Universe evolution), so in the above Eq.© 
9 = 9m + Pde, P = Pm + Pde and we have used the 
unit convention 8ttG/3 — c = 1 throughout. 

As for EOS, taking p — f(p) into Eq.lJSJl, we get gen- 
erally: 



dp f da 

9~W) = ~ 1 ~ 



(4) 



If only we know the concrete form of the function /(p), 
the evolution of scale factor a with the energy density p, 
may be deduced out. 

In order to describing the behavior of dark energy 
with more possibilities, an extension of the Chaplygin 
gas model is considered, and its relatively general EOS 
takes the following form: 



Px 



Ma) 
9x 



-Px 



(5) 



where px is pressure, px (positve) with a > is energy 
density of the dynamically changing part for the dark 
energy component if we also include the cosmological 
constant contribution, and A(a) is a positive function 
depending on scale factor a or its inverse. In some gen- 
eralized Chaplygin gas model the cosmological constant 



can not be obtained at later cosmic times. The total EOS 
of dark energy can be, therefore, divided into two parts — 
the cosmological constant pA and dynamically changing 
part px- That is formerly: 



Pde = 9a + 9x, 
where the interrelations are, 



Pde =Pa+Px 



PA 



-PA, 



PX 



A(a) 
9x 



-9x 



(6) 



(7) 



Also Eq.© is divided into , with the assumptions they 
decouple each other: 



d(p m a 3 ) = -p m d(a 3 ), d(p A a 3 ) = -p A d(a 3 ), 



d(p x a 3 ) = -p x d{a 3 ) 



(8) 



i.e, by assuming that there is no interaction between dif- 
ferent components. In this paper, our main attention is 
put on the 'A-component'. 

When A{a) takes a constant value A and a — 2, it re- 
sults in the following relation, the CG case (see reference 

My 



where B is an integrated constant. We can get a more 
general relationship between density px and scale factor 
a by considering JSJ and © , that is 



a Px,Q 



d(a 3a p a ^ 



3aa 3a - i A{a)da (9) 



where ao is the present value of the scale factor and px,o 
is the present A-term energy density. 

Now entering the essence of this paper, we consider 
the function A(a) with the form A(a) = Aoa~ m (Ao is 
constant, and m > 0), and then the characteristics of our 
model are determined by three undetermined parameters 
(m, a, Aq). A special case with a = 2 results in the 
following relationQ: 



9x 



6 A 



B 



6 — m a m 



Obviously it can easily return to the original CG model 
result by taking m=0 directly. 

Otherwise, in the point of view for intcgrability to 
two cases are obvious as following. We will discuss m — 
3a in the section III Al and m 7^ 3a in the section III Bl 
We have known that present observation data constrains 
the range of the equation of state parameter for dark 
energy as -1.38 < w < -0.82. Hao and Li have 
demonstrated that w = — 1 state is an attractor for the 
Chaplygin gas model and the equation of state of this 
gas could approach this attractor from either w < — 1 or 
w > — 1 sides. So in following sections, we would also 
consider this attractor approaches on the background of 
the extended Chaplygin gas model. 
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FIG. 1: On the pi — a plane, pi changes depending on a. i 
represents M, X, A. (assuming Ao = 1). And The total density 
p = X^i P» i s a l so showed up. 



FIG. 2: On the px — a plane, px changes depending on a. 
With a(t) increasing, px from increases to a positive maxi- 
mum and then decreases to approach a constant near zero. 



The m = 3a case 



In this section, we take m — 3a so Eq.© gives the 
following result: 



Px = ZaA 



,3a 



,3a 



(10) 



We can see that, 



px oc a , a — > ao 

px oc a _3 (lna) i , a — * oo 



It means that around a = ao the dynamically changing 
part dark energy density evolves similar to dust-like mat- 
ter with p ~ a~ 3 . In contrast, it evolves slower as a — > oo 
than that of a ~ ao. 
Figure ^ shows 

from a — 0.5134ao the zero ECG density point, and 
with the increasing of scale factor a, the dynamically 
changing part of dark energy density has reached a 
biggest value (for example: corresponding to a = 2 case, 
{px)max = /Ox(0.6065ao) seeing Figure |2J), and then it 
approaches to a constant near zero. 

As to a — 0, we gain p = —Aop, which can be de- 
scribed directly by conventional EOS p — wp(w is con- 
stant). More discussions for this situation can be found 
in references 

El 13 El, for example. 
Taking Ea. ((TU|) into Fridamman Eq.J2J), we get: 



i^da 



= / dt 



ao (3aA \na/a + a 3 a p Xfi y/^ J to 



Because the uncertainty of a, the above integration is 
not easy to solve. But for a general picture about the 
scale factor, we can simply put a = 1/2. Thus we get 




100 



FIG. 3: The evolution picture of a(t) depends on t, and we 
take a = 1/2. 



the integration result: 
Y 2 Y 3 



\nY + Y 



3A 



2-2! 3-3! 



2a 



3/2 



(11) 

where Y = | In ^ + y / o^x j0 Mo- Eq.JTT} can be directly 
figured as in |3 which explicitly shows scale factor in a 
accelerating state with d > and there is no singularity 
during the process. 

Since no singularity (see reference for the classifica- 
tion of singularity) occurs during the evolution of pxifl), 
we can discuss the w = — 1 crossing with it. So far as we 
know, w = results in the following equations with 
subscript denoting today's value: 

PX I _ 

3aln(a/a ) + l/i) ' 1 A 



W 



PX + PA 
VPX.O + PA 
PX,0 + PA 



(12) 



where i] is defined by rj = Ao/(a m p x ) called Judge for 
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its following property. 

Taking (| 11) [I into above equation, we can get: 



In — I = 



> 



77-1 
3ar] 



provided 




Then we can see that through the evolution of scale fac- 
tor a(t), the value of state parameter w changes contin- 
uously. So that crossing state barrier (w = —1) would 
occur during the procedure of the evolution. It has also 
told that Judge r\ here played an important role of relat- 
ing to the different ranges of the present state parameter 
wo and the scale factor (denoted by a w ) in accordance 
with Phantom Divide(w = — 1) crossing as tabled in lTLAl 
below. 

Table ITTA1 Parameter ranges and w = — 1 



V 


w 


a w (at w = -1) 


= 1 


-1 


= a 


< 1 


> -1 


< a 


> 1 


-1 


> a 



Why can the Judge r] have such important property? 
In the next part, we will work out its physical meaning 
in the more complete background. 

Now, as for eq.©, we can get its equivalent expression 
as: 



h 2 = 



rm w 



k 



H 2 a 2 



PX 

W mfi {l + zf+W X fi{l + zf 
[3a??ln(l + z) + 1]» +W Afi 
+W K0 (l + z) 2 

H(z) 



(13) 



where h(z) with h{z) = £ ' is reduced Hubble pa- 
rameter and z — ao/a — 1 is red-shift, W m ,o = 



W x ,o = ff, W A .o 



3^, and W k<0 



Px,o, Hq, and ao are, respectively, present non-relative 
matter density, dynamically changing part dark energy 
density, Hubble parameter, and scale factor); p C Q = Hq 
is the current critical density. These density parameters 
should satisfy the following relation: 



W m n + Wx.o + W A ,o + W k .o = 1 



(14) 



Previously, we have pointed out that at the stage of 
a — > clq the density px evolves like dust-like matter and 
now we can also see this point from following. Around 
a = ao, we can expand the term [3a7yln(l + z) + l] 1 /" 
around z = 0, and preserve the first order of z and plus 
other terms. Thus, we get: 

h 2 = W m . Q {\ + zf + iW xm {l + z) i + 

+W Afi + W kfi (l + z) 2 (15) 

The equivalent density parameter W m o now is defined 
by: 



W m<0 = W m , + WxA 1 ~ 3v) 



(16) 



which has included the dark energy density contribution 
part. The (z + l) 4 term corresponds to the radiation 
component, which can be neglected in the beginning. 

Considering the deceleration parameter q = —aa/a 2 , 
it can be expressed by the use of reduced Hubble param- 
eter as: 



1 + z d 
2h{z) 2 dz 



Hz) 2 



2(1 + qo)(qo is the present 



q(z) = -1 

Therefore, (dh(zf /dz) z = 
deceleration parameter). Evaluating (dh(z) 2 /dz) 
from Ea. p3|l and solving with respect to ry, we get: 

2q + 2Ha, ~ W m , Q - W x ,o 



(17) 



>1 



3WV 



(18) 



So far, the physical meaning of parameter term rj has 
clearly stood before us. It reflects the combination of den- 
sity parameters and deceleration parameter at present. 
In other words, if we know the values of density pa- 
rameters and deceleration parameter, we would know the 
ranges of wo and a w . In addition, from Eq. l|18|) . we also 
know: 2(q + W A ,o) > W m , + W x ,o- 



B. 



The m/3a case 



In this section, we have m ^ 3a and Eq.@ has the 
following result: 



Px = PxA( 1 



3 a 



-r])x 3a + 



3a 



3a 



3a 



-rjx 



(19) 



where x is defined by a; — ao/a called reduced scale fac- 
tor. We can see that dynamically changing part dark en- 
ergy density evolves depending on two terms of x 3a and 
x m , and the values of pure parameters term — ?y deter- 



mine the effect. Additionally, the equation H19|) also tells 
us that p x decreases while the scale factor a increases. 

Combined with the EOS pde = wPde, the state pa- 
rameter w is expressed as: 



-px + Pa 



w 



Px + Pa 



(20) 



. Taking 1)19(1 into consideration, we can get the following 
relation: 



< 



v> 



(3a — m)/rj — 3a 



w < -1 
provided w = — 1 
w > -1 

Through the evolution of scale factor a, the value of 
x 3a-m would satisfy the above relation, that is, the 



5 



w = — 1 crossing phenomenon would occur. The different 
ranges of the present state parameter wq and the scale 
factor a w corresponding to the state barrier(u) = — 1) 
can be shown from the following table no matter whether 
3a > to or 3a < to. 

Table ITT*Bl Parameter ranges and w = — 1 



V 


wo 


a w (w = -1) 


= 1 


-1 


= ao 


< 1 


-1 


< ao 


> 1 


< -1 


> ao 



Note: it is similar to table ITTa1 but the Judge r] has 
changed. 

Corresponding to Eq. 1|13|) , we can write out the expres- 
sion of the to ^ 3a case as: 

h 2 = T^ m , (l + z) 3 + VI6co(l + ^ 3 [l + 

■"' „((i+zr- 3 --i)]i 



3a — to 



(21) 



It would be interesting to discuss the behavior around 
a = aQ. Thus, we can expand the [1 4- 3 + 
z yn-3a _ 2)]^ arounc j ^ — o, preserve the first order 
of z, and get: 

h 2 = W m , (l + zf- ZWxflVO- + zf + 

+W Afi + W kfl (l + z) 2 ( 22 ) 

where W m ,o is defined by 

W m ,o = W m , + WxM 1 + 3 '?) (23) 

The effective matter component also gets the contribu- 
tion from the dark energy component. 

With the use of Ea. (|17ll . we can get the Judge as: 



2 g0 + 2W A .o ~ W m . - Wx,o 
3W x ,o 



(24) 



It results in:2(g + W Ay o) < W m fl + Wx,o which is quite 
different from to = 3a case, and this difference can be 
used as to distinguishing the two cases. 



III. REMARKS ON A(x) 



series expansion. We substitute A(a) by A(x) for simpli- 
fication, and then expand A{x). 

A ^ = E —^- m 



A(0) + A'(0)x + -A"(0)x 2 + 



-^-A^(0)x n 

m.< 



Taking it into the expression of EOS, the pressure be- 
comes 



\- ^(o) x 

m. ' X 



(25) 



Compared with A{a) — Aoa~ m , the above expression 
is a combination of Eq.JSJ) with different values of to. 
This property may also imply that at the different cos- 
mic stages there is an accordant effective value of m in 
the EOS (refeq:state equation) to describe the dark en- 
ergy fluid. As for far future stage, m takes zero after 
neglecting small quantities. 

Calculating density from Eq.©, we get 



Px 



E 



1 



3a 



i! 3a 



, (o) x 



Often, px and a depend on each other and we can de- 
duce out the interrelationship of them from integration 
formula©. Thus, pressure px can be at last expressed as 
a function with one variable , for example, px(px), but 
not the two variables function px{px,o)- As for EOS 
Px = —^Ir-px, A(x) and a are both used to determine 

Px(px) and at last px{px) has no relation with a. Ac- 
cording to this assumpation, we would try to choose A{x) 
to satisfy this requirment. Through a series of analysis, 
we finally get 



A(x) = B a {l 



\3a-l 



where B is a lower limit for density p x as following. 
Thus, p x becomes 

Px =B{l + Cl a + --- + C^x m + --- + x* a ) a / a = B(l + x) 3 
where CUl = — .I 3 "- 11 ^ and B is the density lower limit 

oa m!(3a — my. J 

defined by 



In the above treatment a,A(a) — Aoa~ m as in the last 
section is just one of the possible forms for A(a). Gener- 
ally, we are not able to obtain A(a) directly from theo- 
retical analysis. However, as to the far future universe in 
which a; (here we also define x — ao/a) will approach to 
zero with redshift negative if we define today's redshift 
value as zero, the behaviour of A(a) is then determined 
by the small quantity x. It is natural to think of Taylor 



lim px = B 

As a result, the EOS 

px = -B*p% (26) 

is irrelevant to a. Then, we will discuss what information 
about universe this EOS gives. 



6 




W m .o, Wkfl,8,nd Wa,,o are the effective density parameters 
and show the X-term behaviors. Noting our prerequisite 
at cosmic far future evolution stage, the reduced Hubble 
parameter can be approximately written as 

h(z) 2 ~ W A>0 

which means the de sitter phase. If our prerequisite is 
not so strict, we can see that for the bigger values of x, 
h(z) 2 ~ W m fl(l + z) 3 that means a dust-like phase. On 
the whole , our universe has experienced a phase transi- 
tion from dust-like stage to a de sitter stage during the 
process of cosmic evolution. 



FIG. 4: It shows us the relationship of a(t) and t. We assume 
parameters ao = 1, to = 0, B = 1. 



A. The expansion of universe 

For the purpose to show the extended Chaplygin gas 
contribution, we can temporally neglect p m and p\ in 
Eq.p|l and assume the universe is flat for simplification. 
Solve Eq.© and get 



1 



x + 1 



V2 = (i-to)Ba (27) 



, 2^1 + x + x + 2 n 
In t= 2 

(3 + 2y/2)x 

where ao and to is the present scale factor and age of 
universe, respectively. This function can easily be shown 
in figure 01 on which the reduced scale factor increases 
accompanying the elapse of cosmic time. And also a > 
implies our universe at the acceleration phase of state. 



B. contributions from dynamical dark energy 

Under the help of Eq.©, reduced Hubble parameter 
h(z) can take the following form: 

h{z) 2 = W mfi {\ + zf + W k ^{\ + zf + 

+Wx,o — (2 + z) 3 + W A . (28) 
Px,o 

As for 'X-term' in above formula, it is easy to deploy 
(2 + z) 3 to get 1 + z, (1 + z) 2 , (1 + z) 3 and constant 
terms. Then combined them with the other terms, we 
can get: 



h(z) 2 = W mfl {l + zf + W kfi (l + z) 2 + 

3B 

+W x ,o (l + z) + W A ,o 

PX.O 



(29) 



where 

~ B ~ 3B 

W mfi = W mfi + W x ,o- — , W k ,o = W kfl + W x ,o- 



PXfi 



W\ = W A .o + Wx,q- 



Px,o 



B 

PX.Q 



C. The w = — 1 crossing problem 

In section two we have discussed the crossing prob- 
lem in general cases, and now we reconsider it with the 
specific situation. Taking EOS into the state parameter 



formula w 



PDE 



, we get: 
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(^PX+PA 
PX + PA 



(30) 



It is obvious that w > — 1, that is to say, there is no cross- 
ing phenomenon in this specific model and this particular 
case belongs to a quintessence-like dark energy. 



IV. CONCLUSION AND DISCUSSION 

In this paper, in order to discuss the properties of 
dark energy with more possibilities, which is used to ex- 
plain currently cosmic accelerating expansion, we have 
extended the Chaplygin gas model by which the Dark 
Fluid concept can be realized with its contributions to ef- 
fective universe compositions. Through a series of math- 
ematical treatment, it shows that there are indeed some 
cases that can realize state parameter w = — 1 crossing. 
Moreover, the physical meanings behind the parameter 
terms like 3A /aQ a p x in Sec lll Al are interpreted. 

Firstly, as for the m — 3a case, we find that px first 
increases a maximum, and then decreases approaching 
to a constant near 0. During the procedure, the pres- 
sure p decreases to — oo. At the neighborhood of a = ag, 
the behavior of the density px is similar to non-relative 
matter with p tx a 3 . With further analysis, we find that 
the state parameter w can cross its traditional barrier 
w = — 1 smoothly without singularity. By the calcula- 
tions of the reduced Hubble parameter h(z) (|15|) . we find 
that the dynamically changing part of dark energy den- 
sity is effectively equivalent to the sum of the matter-like 
part and radiation- like part. After considering the de- 

3,4o ;„ 



celeration parameter q(z), the parameter term 

explicitly expressed as combination of the present den- 
sity parameters W^o (subscript i = m,X,A,k) and the 
present deceleration parameter qo- 
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Secondly, as to the m ^ 3a case, mostly it is similar 
to the m — 3a case. However, different from the former 
case, the density p x is made up of (ao/a) 3 "-term and 
(ao/a) m -term, and it is equivalent to total effects of both 
a non-relative matter and some kinds of dark energy fluid 
with density px = (ao/a) m . Obviously, from Eci. (|19J) px 
always decreases to approach zero with the scale factor 
a(t) increasing. 

Thirdly, after describing the general situations, we also 
discuss a set of specific values of m — 3a — 1 = 0. Then, 
with parameter independent assumpation, concrete cal- 
culations show that the universe is in the state of accel- 
erating expansion with a > 0. At last, we compare our 
model with the typical state parameter w = pde/pde 
obtaining w > — 1, which means that there is never 
w = — 1 crossing phenomenon in this case. By discussing 
the reduced Hubble parameter h(z), we know that the 
dark energy fluid is evolving equivalently to a fluid with 
px = 3£(i + A). 

To sum up, we have discussed the extended Chaply- 



gin gas model in detail, with the hope to show more 
possible properties for dark energy and to realize the 
concept of Dark Fluid for a suggestive describing the 
cosmic dark components. To constrain the EOS char- 
acterizing parameters, work will be done to maximize 
the following likelihood function(see reference [23j):L oc 
exp[— x 2 (-P)/2]. And also in this paper, interactions be- 
tween the different components of cosmology are not con- 
sidered. We will publish the contents elsewhere soonp^|. 
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